'''
Created on 12.09.2012

@author: Julia
'''

from SpinArt import *
from Physics import *

diamond=LatticeType()
diamond.load_latticetype(0)

nv=SpinType()
nv.load_art(0, "Spins.txt")
c13=SpinType()
c13.load_art(1, "Spins.txt")

lat=OneLattice(diamond)
centralspin=OneSpin(nv)
lat.add_spin(centralspin)
lat.create_spin(c13, 100, 1.5, 1)
lat.subDict["13C"].set_concentration(0.3)

n15=SpinType()
n15.load_art(3, "Spins.txt")

n14=SpinType()
n14.load_art(2, "Spins.txt")


n14.zfs=-5.04
n14.gyro=760*0.4038

nv.gyro=-2.0023*9.2740154e-24/sconst.h/1e4
nv.zfs=2871.5

#print "c13 concentration:", lat.subDict["13C"].conc
#print "Anzahl c13:", lat.subDict["13C"].n_spins
#print "Anzahl max:", lat.subDict["13C"].n_max
#print "Anzahl NVs:", lat.subDict["NV"].n_spins

